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ABSTRACT 

We use the statistics of strong gravitational lenses to investigate whether mass 
profiles with a flat density core are supported. The probability for lensing by ha- 
los modeled by a nonsingular truncated isothermal sphere (NTIS) with image 
separations greater than a certain value (ranging from 0" to 10") is calculated. 
NTIS is an analytical model for the postcoUapse equilibrium structure of viri- 
alized objects derived by Shapiro, Iliev, & Raga. This profile has a soft core 
and matches quite well with the mass profiles of dark matter-dominated dwarf 
galaxies deduced from their observed rotation curves. It also agrees well with the 
NFW (Navarro-Frenk- White) profile at all radii outside of a few NTIS core radii. 
Unfortunately, comparing the results with those for singular lensing halos (NFW 
and SIS-I-NFW) and strong lensing observations, the probabilities for lensing by 
NTIS halos are far too low. As this result is valid for any other nonsingular 
density profiles (with a large core radius), we conclude that nonsingular density 
profiles (with a large core radius) for CDM halos are ruled out by statistics of 
strong gravitational lenses. 

Subject headings: cosmology: theory - cosmology: observations - gravitational 
lensing 

1. Introduction 

The concordant, currently dark energy-dominated, spatially fiat cold dark matter (CDM) 
cosmology is so successful that we are now said to be in an era of "precision cosmology" (Pee- 
bles 2002; Ostriker & Souradeep 2004). In this concordant cosmology, dark energy (DE), 
mainly introduced to close the universe and to explain the accelerating expansion of the 
universe, attracts much attention in the realms of astrophysics and theoretical physics. The 
simplest candidate for DE is the cosmological constant, and its more general form, known 
as quintessence, is a cosmic scalar field minimally coupled with the usual matter (Caldwell, 
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Dave, & Steinhardt 1998; Ma et al. 1999; Peebles & Ratra 2003). The generalized Chap- 
lygin gas, as a unification of dark matter and dark energy, was recently proposed and can 
be constrained by observations (Padmanabhan & Choudhury 2002; Zhu 2004). In other 
cosmological models, DE is replaced by certain possible mechanisms, such as brane world 
cosmologies (Randall & Sundrum 1999a,b) and the Cardassian expansion model (Frcese & 
Lewis 2002; Zhu & Fujimoto 2002, 2003, 2004). Despite its success, the standard CDM 
theory of cosmic structure formation has several problems, which exist mostly in the small- 
scale regime. For example, the observed rotation curves of dark-matter-dominated dwarf 
and low surface brightness (LSB) disk galaxies tend to favor mass profiles with a fiat density 
core (e.g., Salucci & Burkert 2000; Gentile et al. 2004), unlike the singular profiles of the 
CDM A^-body simulations (e.g., Navarro, Frenk & White 1997, NFW profiles; Moore et al. 
1999; Jing 2000; Jing & Suto 2002) and that favored by baryon cooling models (i.e., singular 
isothermal sphere, SIS profile). While there are debates on whether the observed data were 
resolved well enough to indicate a soft core (van den Bosch & Swaters 2001; Marchesini et 
al. 2002), quite recent A^-body simulations of CDM with higher and higher force and mass 
resolution still favor cuspy halo profiles (Diemand, Moore & Stadel 2004; Fukushige, Kawai 
& Makino 2004; Navarro et al. 2004; Tasitsiomi et al. 2004; Wambsganss, Bode & Ostriker 
2004). 

Recently, an analytical model was presented for the post-collapse equilibrium structure 
of virialized objects that condense out of a low-density cosmological background universe, 
with or without a cosmological constant (Shapiro, Ihev & Raga 1999; Iliev & Shapiro 2001). 
The model is based on the assumption that cosmological halos form from the collapse and 
virialization of 'top-hat' density perturbations, and are spherical, isotropic and isothermal. 
According to the authors, this predicts a unique, non-singular, truncated isothermal sphere 
(NTIS) and provides a simple physical clue about the existence of soft cores in halos of 
cosmological origin. This NTIS model is claimed to be in good agreement with observations of 
the internal structure of dark matter dominated halos on scales ranging from dwarf galaxies to 
X-ray clusters. In particular, it matches quite well the mass profiles of dark matter dominated 
dwarf galaxies deduced from their observed rotation curves (Shapiro et al. 2004). Quite 
recently, the NTIS model was revisited by using the self-interacting dark matter (SIDM) 
hypothesis (Ahn & Shapiro 2004). There are other efforts for analytical derivations of the 
density profile; e.g., Hansen (2004) derived the bound on the central density slope of -1 
analytically (as found numerically by NFW). This is done by a simple solution to the Jeans 
equation, which is valid under the assumption that both the central density profile and the 
phase-space-like density are exact power laws. However, this work did not clearly give a 
density core. 

To investigate whether there is a soft core at the center of each CDM halo, we use another 
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independent and robust tool, gravitational lensing, which has been widely used to detect the 
mass distributions in the Universe (Turner, Ostriker, & Gott 1984; Schneider, Ehlers, & 
Falco 1992; Wamsganss et al. 1995; Wu 1996; Bartelmann 1996; Wu 2000; Bartelmann & 
Schneider 2001; Chen 2001; Keeton 2001; Keeton & Madau 2001; Xue & Wu 2001; Keeton 
2002, 2003; Pen et al. 2003; Schneider 2003; Keeton & Zabludoff 2004; Wu 2004; Zhang 
et al. 2004) and the dark energy density and its equation of state (Fukugita, Futamase, & 
Kasai 1990; Fukugita, & Turner 1991; Turner 1990; Krauss & White 1992; Maoz & Rix 1993; 
Kochanek 1995, 1996; Falco, Kochanek, & Muiioz 1998; Cooray & Huterer 1999; Waga & 
Miceh 1999; Sarbu, Rusin, & Ma 2001; Dev, Jain & Alcaniz 2004). Recently, motivated by 
the current largest homogeneous sample of lensed quasars, coming from the radio Cosmic 
Lens All-Sky Survcy(CLASS; Myers et al. 2003; Browne et al. 2003), an extension of the 
earlier Jodrell Bank/Very Large Array Astromctric Survey (JVAS; Patnaik ct al. 1992; 
King et al. 1999), much work has been devoted to statistics of strong lensing to constrain 
the cosmological parameters and mass distribution (Chae et al. 2002; Li & Ostriker 2002; 
Oguri et al. 2002; Chae 2003; Li & Ostriker 2003; Oguri 2003; Oguri, Suto k Turner 2003; 
Oguri & Keeton 2004; Wang 2004; Mitchell et al. 2005; Sereno 2005; Zhang et al. 2005). 
Statistics of strong gravitational lensing by halos with a soft core was studied quite early 
(Hinshaw & Krauss 1987; Chiba & Yoshii 1999; Cheng & Krauss 2000), but no observational 
sample suitable for analysis could be used in their work; in particular, they did not use their 
cored density profile to fit the observed rotation curves of the dwarf and LSB disk galaxies. 
In strong gravitational lensing of quasars, the separation between multiple images is the 
most important observable. In some instances, asymmetry and shear of lens halos affect the 
properties of images considerably, but statistically, image separations are mainly determined 
by the potential wells of the host halos characterized by the slopes of their density profiles. 
Thus, the probabilities for lensing by NTIS halos with image separations greater than 
provide us with an independent and powerful probe of the existence of soft cores of CDM 
halos. Our calculations are performed in a concordant, spatially flat ACDM cosmology 
favored by the Wilkinson Microwave Anisotropy Probe ( WMAP; Bennett et al. 2003) plus 
large-scale structure (LSS) data from the Sloan Digital Sky Survey (SDSS; Tegmark et al. 
2004a,b), with matter density parameter = 0.3, galaxy fluctuation amplitude as = 0.9, 
and Hubble parameter h = 0.72. For comparison, we also give the results when NTIS is 
replaced by SIS+NFW and NFW proflles. 



2. Lensing Probabilities 

In what follows, we use a reduced mass of a halo deflned as M15 = M/(lO^^h^-'^M0). 
The NTIS model is based on the assumption that cosmological halos form from the collapse 
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and virialization of top-hat density perturbations and are spherical, isotropic, and isothermal 
(Shapiro, Iliev & Raga 1999; Iliev & Shapiro 2001). The well-fitted density profile is given 
by (Shapiro, Iliev & Raga 1999; Ihev & Shapiro 2001; Shapiro et al. 2004) 
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where A — 21.38, B — 19.81, a — 3.01, and b — 3.28. The central value of the density profile 
is Po = 1-8 X 10^pc(-2coii), where Pc(^coii) is the critical density of the Universe at the epoch 
of halo collapse with redshift 2;coii- The small core radius depends on mass M and 2;coii and 
is given by 

ro = 0.115(M/po)'/' = 6.73 x lO-^M^f [1^^(1 + ^eou)' + n^]-'/^h-'Mpc, (2) 

where we have used Pc{z) = pc(0)[Jlm(l + z)^ + ^a] (with fl^ = 0.3 and Qa = 0.7 in 
later calculations) and pc(0) = 1.88 x IQ-'^^h'^gcm.-^ = 2.777 x W^H'^MqUpc'^. It should 
be pointed out that the original formula for tq, equation (88) in Shapiro et al. (2004), is 
To = 1.51 X 10^'^(M/po)"^^^ (we write M200 as M, which is discussed in section 3). This is 
a clerical error, but we have checked that their subsequent results are not affected by this 
wrong formula. 

The gravitational lens equation is 77 = Dg^/Di^ — Di^^a, where 7] and ^ arc the physical 
positions of a source in the source plane and an image in the image plane, respectively, a is the 
deflection angle, and Dl, D^, and Dls are the angular diameter distances from observer to 
lens, observer to source, and lens to source, respectively. By defining dimensionless positions 
y — Di^Tj/DsTo and x — ^/ro, and dimensionless angle a — Di^Di^^a/ D^tq, the lens equation 
is then (Shapiro et al. 2004) 

y^x- , ^'^^'^l , f AVa2 + x2 - bVWT^ -Aa + Bb) , (3) 
(^40 — Ba) \ / 

where Kc is the central convergence: 

s(e ^ 0) WoM B\ 
i^c = = X ' 

^crit ^crit \ ^ / 

where E^it = c^Ds/^nGDi,Di^s is the critical surface density. 

It is well known that generally, for any spherically symmetric density profiles of lensing 
halos, multiple images can be produced only if the central convergence is greater than unity 
(Schneider, Ehlers, & Falco 1992). When < 1, only one image is produced. Note that 
even ii Kc> 1 is satisfied, multiple images can occur only when the source is located within 
Ucr — yixci), where x^r is determined from the lensing equation (eq. [3]) with dy/dx — 
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for X < (this is similar to lensing by NFW halos). For a singular density profile such as 
the SIS and NFW profiles, the central value is divergent, so k > 1 is always satisfied, and 
multiple images can be produced for any given mass. For density profiles with a finite soft 
core such as the NTIS profile, however, the condition k, > 1 requires that only halos with 
mass greater than a certain value (determined by Kc = 1) can produce multiple images. This 
is clearly shown in Figure 1, where three curves for Kc = 1.1, 1.05, and 1.0 are plotted, and 
when Kc — 1.0, only one image is produced. In lensing statistics, this requirement will limit 
the populations of lensing halos to quite a small fraction. Such a conclusion is valid for any 
lensing halos with a finite soft core, which is discussed in detail later. 

When quasars at redshift Zg are lenscd by foreground CDM halos of galaxies and clusters 
of galaxies, the lensing probability for image separations larger than is (Turner, Ostriker, 
& Gott 1984; Schneider, Ehlers, & Falco 1992) 

P{> A0) = / V{z,)dzs / -^^dz / n{M, z)a{M, z)B{M, z)dM, (5) 
J Jo dz Jq 

where V{zs) is the redshift distribution for quasars approximated by a Gaussian model with 
a mean of 1.27 and a dispersion of 0.95 (Helbig et al. 1999; Marlow et al. 2000; Chae et al. 
2002; Myers et al. 2003), Df^{z) is the proper distance from the observer to the lens located at 
redshift z, n{M, z) is the physical number density of virializcd dark halos of masses between 
M and M+dM, and B{M, z) is the magnification bias. The physical number density n(M, z) 
is related to the comoving number density n(M, z) by n(M, z) = n{M, z){l + z)^] the latter 
is originally given by Press & Schechter (1974), and the extended version is n{M,z)dM — 
{po/M)f{M, z)dM, where po is the current mean mass density of the universe and 

f{M,z) = (0.315/M)(dlnA,/dlnM)exp(-|ln(A,/1.68) + 0.61p-^) (6) 

is the mass function for which we use the expression given by Jenkins et al. (2001). In 
this expression, A^ = 5c{z)/A{M), in which Sc{z) is the overdensity threshold for spherical 
collapse at redshift z and A(M) is the rms of the present variance of the fluctuations in a 
sphere containing a mean mass M. The overdensity threshold is given by Sc{z) = 1.68/ D[z) 
for the ACDM cosmology (Navarro, Frenk, & White 1997), where D{z) = g[Q{z)]/[g{Q^){l + 
z)] is the linear growth function of the density perturbation (Carroll, Press, & Turner 1992), 
in which g{x) = 0.5a;(l/70 + 209a;/140 - a;Vl40 + x'^^^)-^ and n{z) = ^^(1 + zf/[l - 
+ ^^m(l + z)^]- When we calculate the variance of the fluctuations A^(M), we use 
the fitting formulae for the CDM power spectrum P{k) — AkT'^{k) given by Eisenstein 
& Hu (1999), where A is the amplitude normahzed to (jg = A(rM = 8/i~^Mpc), given by 
observations. The cross section for lensing is a{M, z) = 7iy^^r'^'d{M — Mmin), where i?(a;) is a 
step function, and Mmin is the minimum mass of halos above which lenses can produce images 
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with separations greater than A^. The minimum mass Mmin can be derived directly from 
the relationship between the mass of a lens halo and the corresponding image separation, 
as follows. It is obvious from Figure 1 that the separation between the outer two images is 
almost independent of the source position y. Thus, similar to the analysis for NFW profiles 
(Li & Ostriker 2002), the image separation Ax{y) produced by a NTIS halo for a source at y 
can be well approximated by Ax{0) = 2xo, where xo is the Einstein radius determined from 
the lens equation with y = 0. The angular image separation is 

Axro 2xoro 

Then, from this equation and equation (2) we have 

M^^ = 1.26 X 10^2 1^1^ (^^^ (J^^ ' [^^(1 + + f]^]h-iM0, (8) 

where c/ Hq — 2997.9/i~-'^Mpc is the present Hubble radius. The magnification bias B{M, z) 
is calculated numerically with 
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(M, ^) = (2/y,,) / dyy\ii{y)Y-' (9) 
Jo 



given by Oguri et al. (2002), where ^{y) is the total magnification of the two outer images 
for a source at y; this can also be computed numerically. 

The numerical results of Eq. (5) for NTIS lens halos are plotted in figure 2 {thin solid 
line). For comparison, the survey results of JVAS/CLASS and the predicted probability for 
lensing by SIS + NFW and NFW profiles are also shown. A subset of 8958 sources from the 
combined JVAS/CLASS survey form a well-defined statistical sample containing 13 multiply 
imaged sources (lens systems) suitable for analysis of the lens statistics (Myers et al. 2003; 
Browne et al. 2003; Patnaik et al. 1992; King et al. 1999). The observed lensing probabilities 
can be easily calculated (Chen 2003b, 2004a) by Pobs(> A^) = A^(> Ae)/8958, where A^(> 
A^) is the number of lenses with separation greater than A^ in 13 lenses. The observational 
probability Pobs(> A.9) is plotted as a histogram in Figure 2. In the two-population model 
SIS + NFW, the galaxy-size and the cluster-size lens halos are approximated by SIS and 
NFW profiles, respectively (Sarbu, Rusin, & Ma 2001; Li & Ostriker 2002; Chen 2003a,b, 
2004a,b; Zhang 2004). In the one-population model NFW, lens halos with different sizes arc 
all approximated by NFW profiles (Li & Ostriker 2002); this is similar to the NTIS model. 
The theoretically predicted lensing probabilities shown in Figure 2 are calculated separately 
for three cases (NTIS, SIS-^NFW and NFW) according to equation (5) . The differences in 
lensing probability distributions for the three cases arise from their different values for the 
lensing cross section (t(M, z) and magnification bias B{M, z), since these two quantities are 
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determined uniquely from the corresponding density profile . Here we have recalculated the 
Icnsing probabilities for SIS + NFW and NFW profiles according to equation (5). Since the 
density profiles, lensing equations, and lensing cross sections for SIS and NFW profiles have 
been discussed many times in the literature, here we only give the final results of lensing 
probabilities for these two models (for details about SIS+NFW and NFW profiles, see Li & 
Ostriker 2002; Chen 2004a, and references therein). 



3. Discussion and Conclusions 

One can see clearly from Figure 2 that probability for lensing by NTIS halos with 
an image separation greater than is far too low to match the observational results of 
CLASS/ JVAS; it is even much less than that for NFW lenses. We thus conclude that, at 
least, NTIS as a model to approximate density profiles of dark matter halos is ruled out by 
statistical strong lensing. Within the framework of statistics of strong lensing as displayed 
in the literature, no mechanism can save this model. 

In fact the above conclusion is general for any mass profile with a fiat soft core charac- 
terized by the core density pcore and core radius rcore (these two parameters are determined 
by rotation curves of dark matter dominated dwarf and LSB disk galaxies). To see this, we 
revisit another such density profile, an isothermal sphere with a soft core (cored isothermal 
sphere, CIS): p{r) = al/2TTG{r'^+rl^,^) (Hinshaw & Krauss 1987; Chiba & Yoshii 1999; Cheng 
& Krauss 2000), where Uy is the one-dimensional velocity dispersion. The lens equation is 

y = x- 2nf^{^x^ + 1 - l)/x, (10) 

where y and x are defined in the same way as in equation (3) (here x — C/T^core) and 

S = a2/2G'reoreScnt (11) 

is the central convergence. Similarly, multiple images can be produced if and only if > 1. 
The corresponding lens equation is plotted in Figure 3 for three different values of k^^; the 
curves are quite similar to the NTIS model, even though their density profiles seem quite 
different. 

It is not difficult to understand the extremely low value of the probability for lensing by 
dark halos with a flat soft core. In fact, in our previous calculations, for the NTIS profile, 
Kc can be written in the form 

= 475.7^^^^M'£ = 3.65M,f [0^(1 + zf + ^k?"^. (12) 
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where — Di^Di^s/Ds ■ We have defined, as usual, the mass of a dark matter halo to be 
M — M200 = 4:71 f^^°° p{r)r'^dr (r2oo has its usual definition) and assume that the redshift ^cou 
by which the dark halo collapsed is equal to the lens location redshift and that such an 
assumption will not affect our results (Shapiro et al. 2004). For a source at Zs = 3.0 and the 

1 /3 

lens at = 0.5, equation (12) gives Kc = l.OSM^g . In this case, Kc = 1 implies M ~ M15, 
and this means that for this typical lens system, multiple images can be produced only 
when the lens mass is higher than lO^^h^^Mg, which is the typical size of galaxy clusters. 
By setting Kc — 1 in equation (12), the minimum lens mass for producing multiple images 
is determined as a function of both and Zg. In Figure 4, we plot the mass M{kc — 1) 
(measured in M15) as a function of Zs for three given values of zi,: 0.1, 0.5, and 1.0. For a given 
^L, the minimum mass decreases with increasing Zg. Since the mean value of the redshift of 
the sources is (^s) = 1.27 in our calculations, the actual mean value of the mass for producing 
multiple images (for lenses at Zi^ = 0.5) is larger than lO^^h~^M0. It is also obvious in Figure 
4 that in most cases, multiple images can be produced only when the halo mass is larger 
than lO^'^h^^MQ. From equation (5), we know that the lensing probability is determined 
by the number density of dark halos and the cross section. As pointed out by Shapiro et 
al. (2004), an NFW profile would produce more multiple-image lenses than an NTIS one at 
relatively lower masses, and this trend is reversed at higher masses. Statistically, however, it 
is the requirement of /tc > 1 and thus the existence of a large core radius that strongly limits 
the populations of lensing halos [number density n{M, z)] to quite a small fraction. Thus, 
the extremely low value of the probability for lensing by NTIS halos arises from the quite 
low value of the corresponding number density of galaxy clusters. Furthermore, the lensing 
cross section a ~ 7/^^,, and from Figure 1 and Figure 3 we see that i/cr is quite sensitive to 
— Kc{> 1) — 1 rather than Kc itself. Therefore, a shght change in Kc ~ M^^^ will result 
in a large change both in and xq (the Einstein radius), and since xq ~ A^, this explains 
the quite fiat curve for the NTIS lensing probabihty in Figure 2 {thin solid line). Namely, 
the insensitivity of NTIS lensing probabihty to reflects the fact that within the image 
separations of a few arcseconds, the lens mass and thus the corresponding number density 
of dark halos around that mass change only a little. Similarly, for a CIS profile, after solving 
the equations 

M = 8007rpcr2oo/3 (13) 

- 47r / drr'al/27rG{r' + rlJ (14) 
Jo 

~ 2(7^(r200 - Trrcore/S)/^ (15) 

for r2oo and cr^, can be related to dark halo mass M. The above analysis for the NTIS 
profile then is also true for the CIS model, if core radius rcore is large enough to be able to 
fit the observed rotation curves of the dwarf and LSB disk galaxies. 



-9- 



Note that our conclusion that a CIS model as a cored mass profile is ruled out by 

statistical strong gravitational lensing seems inconsistent with previous works (Hinshaw & 
Krauss 1987; Chiba & Yoshii 1999; Cheng & Krauss 2000), since these authors used this 
model to constrain cosmological parameters. A simple analysis, however, shows us that 
there is no discrepancy. The central convergence can be expressed in terms of the mass 
and the core radius of a lens halo. For NTIS, from equation (2) and equation (4), we 
have Kc oc M/tq] for CIS, from equation (11), equation (13) and equation (15) we have 
(approximately) oc M^/^/rcore- A larger Tcore needs a larger M to ensure > 1, or, 
for a fixed value of mass M, an appropriate smaller value of rcore would ensure > 1. 
Therefore, if Tcore is not determined by currently observed rotation curves of the dwarf and 
LSB disk galaxies, but rather is adjustable (Hinshaw & Krauss 1987; Chiba & Yoshii 1999; 
Cheng & Krauss 2000), then no discrepancy appears. A more realistic mass profile should 
not be divergent at the center (cusp), that is, it should have a fiat core, but the core radius 
Tcore should be small enough to ensure that the lensing probability matches the observations 
of CLASS/JVAS. 

I thank the anonymous referee for quite useful suggestions. This work was supported 
by the National Natural Science Foundation of China under grant 10233040. 
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Fig. 1. — Lens equation for NTIS halos, plotted according to equation (3). Coordinates y 
and X are dimensionless positions of the source and the image, respectively. Curves for three 
values of Kc —1.1 {solid curve), 1.05 {dotted curve) and 1.0 {dashed curve) are plotted. When 
Kc — 1, only one image is created. 



- 15 - 




I I \ I I I I I I I I 

2 2 4 6 8 10 

A6(arc seconds) 

Fig. 2. — Predicted probability for an image separation angle greater than for lensing by 
NTIS halos {thin solid line). For comparison, the survey results of JVAS/CLASS (thick his- 
togram), the predicted probability for lensing by SIS+NFW {dotted line) and NFW {dashed 
line) profiles are also shown. 
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Fig. 3. — Lens equation for CIS halos, plotted according to equation (10). Coordinates y 
and X are dimensionless positions of the source and the image, respectively. Curves for three 
values of =1.1 {solid curve), 1.05 {dotted curve) and 1.0 {dashed curve) are plotted. When 
Kc — 1-1 only one image is created. 
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Fig. 4. — Lens mass M{kc = 1) as a function of ^s, determined by setting = 1 in equation 
(12) for given values of zi^. Curves for three typical values of Zi, —0.1 {solid curve), 0.5 
{dashed curve) and 1.0 {dotted curve) are plotted. 



